Brittle/quasi-brittle transition in dynamic fracture: An energetic signature 
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Dynamic fracture experiments were performed in PMMA over a wide range of velocities and reveal 
that the fracture energy exhibits an abrupt 3-folds increase from its value at crack initiation at a 
well-defined critical velocity, below the one associated to the onset of micro-branching instability. 
This transition is associated with the appearance of conies patterns on fracture surfaces that, in 
many materials, are the signature of damage spreading through the nucleation and growth of micro- 
cracks. A simple model allows to relate both the energetic and fractographic measurements. These 
results suggest that dynamic fracture at low velocities in amorphous materials is controlled by the 
brittle/quasi-brittle transition studied here. 

PACS numbers: 46.50.+a, 62.20.M-, 61.43.-j 



Dynamic fracture drives catastrophic material failures. 
Over the last century, a coherent theoretical framework, 
the so-called Linear Elastic Fracture Mechanics (LEFM) 
has developed and provides a quantitative description of 
the motion of a single smooth crack in a linear elastic 
material [ij. LEFM assumes that all the mechanical en- 
ergy released during fracturing is dissipated at the crack 
tip. Defining the fracture energy F as the energy needed 
to create two crack surfaces of a unit area, the instan- 
taneous crack growth velocity v is then selected by the 
balance between the energy flux and the dissipation rate 
Tv. This yields Q: 



{1-v/cr)K\c)/E, 



(1) 



where cr and E are the Rayleigh wave speed and the 
Young modulus of the material, respectively, and K{c) is 
the Stress Intensity Factor (SIF) for a quasi-static crack 
of length c. K depends only on the applied loading and 
specimen geometry, and characterizes entirely the stress 
field in the vicinity of the crack front. 

Equation ([!]) describes quantitatively the experimental 
results for dynamic brittle fracture at slow crack veloc- 
ities However, large discrepancies are observed in 
brittle amorphous materials at high velocities In 
particular (i) the measured maximum crack speeds lie in 
the range 0.5 — 0.6ci?, i.e. far smaller than the limiting 
speed Cr predicted by Eq. ([1]) and (ii) fracture surfaces 
become rough at high velocities (see [s!, I3| for reviews). 
It has been argued [7] that experiments start to depart 
from theory above a critical ^Acr associated to the 
onset of micro-branching instabilities [8]: for v > vi) the 
crack motion becomes a multi-cracks state. This trans- 
lates into (i) a dramatic increase of the fracture energy 
F at due to the increasing number of micro-branches 
propagating simultaneously and (ii) a non-univocal rela- 
tion between F and v . The micro-branching instability 
hence yielded many recent theoretical efforts [9]. How- 



ever, a number of puzzling observations remain at smaller 
velocities. In particular, even for velocities much lower 
than vi)^ (i) the measured dynamic fracture energy is gen- 
erally much higher than that at crack initiation 
[12I and (ii) fracture surfaces roughen over length scales 
much larger than the microstructure scale ("mist" pat- 
terns) [13], the origin of which remains debated fl3, flBt. 

In this Letter, we report dynamic fracture experiments 
in polymethylmethacrylate (PMMA), the archetype of 
brittle amorphous materials, designed to unravel the pri- 
mary cause of these last discrepancies. We show that 
dynamic fracture energy exhibits an abrupt 3-folds in- 
crease from its value at crack initiation at a well-defined 
critical velocity Va well below ^5. This increase coincides 
with the onset of damage spreading through the nucle- 
ation and growth of micro-cracks, the signature of which 
is the presence of conic patterns on post-mortem frac- 
ture surfaces. A simple model for this nominally brittle 
to quasi-brittle transition is shown to reproduce both the 
energetic and fractographic measurements. 

Dynamic cracks are driven in PMMA with measured 
Young modulus and Poisson ratio of E = 2.8±0.2 GPa 
and u = 0.36, which yields cr = 880 ± 30 m.s~^. Its 
fracture energy at the onset of crack propagation was de- 
termined to be K^/E = 0.42 ± 0.07 kJ.m-^, with Kc 
being the material toughness. Specimen are prepared 
from 140 x 125 x 15 mm^ parallelepipeds in the x (prop- 
agation), y (loading) and z (thickness) directions by cut- 
ting a 25 X 25 mm^ rectangle from the middle of one 
of the 125 X 15 mm^ edges and then cutting a 10 mm 
groove deeper into the specimen (Fig. 1, bottom inset). 
Two steel jaws equipped with rollers are placed on both 
sides of the cut-out rectangle and a steel wedge (semi- 
angle 15°) is pushed between them at constant velocity 
38 /im.s~^ up to crack initiation. In this so-called wedge 
splitting geometry, the SIF K decreases with the crack 
length c. To increase its value at crack initiation, and 
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FIG. 1: Measured crack velocity v as a function of crack 
length c in a typical experiment {Uo = 2.6 J). The verti- 
cal lines are error bars. Top inset: Calculated quasi-static 
SIF K as a function of c. Bottom inset: Schematics of the 
Wedge- Splitting test. 



therefore the initial crack velocity, a circular hole with a 
radius ranging between 2 and 8 mm is drilled at the tip of 
the groove to tune the stored mechanical energy Uq. Dy- 
namic crack growth with instantaneous velocities ranging 
from 75 m.s~^ to 500 m.s~^ and stable trajectories are 
obtained. The location c{t) of the crack front is measured 
during each experiment (40 /im and 0.1 /is resolutions) 
using a modified version of the potential drop technique: 
A series of 90 parallel conductive lines (2.4 nm-thick Cr 
layer covered with 23 nm-thick Au layer), 500 /im-wide 
with an x-period of 1 mm are deposited on one of the x-y 
surfaces of the specimen, connected in parallel and ali- 
mented with a voltage source. As the crack propagates, 
the conductive lines are cut at successive times, these 
events being detected with an oscilloscope. The instan- 
taneous crack velocity v{c) is computed from c(t), and 
the instantaneous SIF K{c) is calculated using 2D finite 
element calculations (software Castem 2007) on the exact 
experimental geometry, assuming plane stress conditions 
and a constant wedge position as boundary condition. 

Values for the fracture energy F are obtained directly 
from Eq. ([!]) by combining the v measurements and the 
K calculations. Typical v{c) and K{c) curves are shown 
in Fig. [H The variations of F with v (Fig. [2]) are found 
to be the same in various experiments performed with 
various stored mechanical energy Uq > 2.0 J at crack ini- 
tiation. This curve provides evidence for three regimes, 
separated by two critical velocities. For slow crack ve- 
locities, F remains of the order of K'^/E as expected 
in LEFM. Then, as v reaches the first critical velocity 
Va — 165m. = O.lQci?, F increases abruptly to a 
value about 3 times larger than K^jE. Beyond Va^ F 
increases slowly with v up to the second critical veloc- 
ity, v\) = 0.36cr 317 m.s~^ 0, above which F di- 
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FIG. 2: (color online). Fracture energy F as a function of 
crack velocity v for five different experiments with different 
stored mechanical energies Uo at crack initiation: 2.0 (□), 2.6 
(o), 2.9 (0), 3.8 (+) and 4.2 J (x). The two vertical dashed 
lines correspond to Va and Vb . The two horizontal dashed lines 
indicate the confidence interval for the measured fracture en- 
ergy Kc/E at crack initiation. Thick red line: model predic- 
tion. Inset: F as a function of K'^/E (see model) for the same 
experiments. A crossover between two linear regimes (linear 
fits in black lines) occurs at {K^/E ^ KljE kJ.m~^ ; 

F = Fa ^ 1.34 kJ.m-2). 



verges again with v. This second increase corresponds to 
the onset of the micro-branching instability, widely dis- 
cussed in the literature 0, 111, whereas the first one, at 
is reported here for the first time. The high slope of 
V{v) around Va provides a direct interpretation for the 
repeated observation of cracks that span a large range of 
F but propagate at a nearly constant velocity of about 
0.2ci^ (see e.g. refs. [ii,[i3|). 

The post-mortem fracture surfaces shed light on the 
nature of the transition at = on the curve r(v). 
Fig. [3] shows the surface morphology for increasing crack 
velocity. For v < Va, the fracture surfaces remain smooth 
at the optical scale (Fig. [31(a), top). Above Va conic 
marks are observed (Figs. [31(b) and [31(c), top). They do 
not leave any visible print on the sides of the specimens 
(Fig. [31(b), bottom), contrary to the micro-branches that 
develop for v > vjj (Fig. [Sjc), bottom). 

Similar conic marks were reported in the fracture 
of many other amorphous brittle materials (see 0, 
and references therein), including polymer glasses, silica 
glasses and poly crystals. Their formation is thought to 
arise from inherent toughness fluctuations at the micro- 
structure scale due to material heterogeneities randomly 
distributed in the material (itI . [i8| . The enhanced stress 
field in the vicinity of the main crack front activates some 
of the low toughness zones and triggers the initiation of 
secondary penny-shaped micro-cracks ahead of the crack 
front. Each micro-crack grows radially under the stress 
associated with the main crack along a plane different 
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FIG. 3: Microscope images (xlO) taken at (a) v — 120 ± 
20 m.s-\ K^/E = 1 (b) v = 260dz30 m.s-\ K^/E = 

2 kJ.m-2 (c) = 650 zb 100 m-s"^ {K'^/E = 7 kJ.m-^). Top 
line : fracture surfaces (0.5 x 0.7 mm^ field of view). Bottom 
line : sample sides (0.25 x 0.7 mm^ field of view). Crack 
propagation is from left to right. 



from it. When two cracks intersect in space and time, 
the ligament separating them breaks up, leaving a visi- 
ble conic marking on the post-mortem fracture surface. 

Figure [4] shows the surface density of conic marks p as 
a function of crack velocity v. Below Va, no conic mark 
is observed up to x50 magnification, consistently with 
(l9| . Above Va, p increases almost linearly with v — Va- 
The exact correspondence between the critical velocity 
Va at which F exhibits an abrupt increase and the veloc- 
ity at which the first conic marks appear on the fracture 
surfaces strongly suggests that both phenomena are as- 
sociated with the same transition. The nucleation and 
growth of micro-cracks can therefore be identified as the 
new fracture mechanism that starts at Va- This damage 
process is generic in brittle materials and is relevant for 
an even wider range of materials than those that exhibit 
conic marks, e.g. granite [20]. 
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FIG. 4: (color online). Surface density p of conic marks as a 
function of crack velocity for all experiments shown in Fig. [2l 
Inset: p as a function of K'^/E (linear fit in black line). 



We now present a simple model reproducing the r{v) 
curve between and vi,. We assume that linear elas- 
ticity fails in the material when the local stress reaches 



a yield stress ay- It defines a fracture process zone 
(FPZ) around the crack tip, the size of which is given 
by Rciv) = v)/a(jY where a is a dimensionless con- 

stant [l5| and Kd is the dynamic SIF. We consider that all 
the dissipative phenomena (plastic deformations, crazing 
or cavitation for instance) occur in the FPZ, with a vo- 
lumic dissipated energy e. The material is then assumed 
to contain a volume density ps of discrete "source-sinks" 
(SS, see e.g. [HI for previous uses of this concept). Each 
SS is assumed to activate into a micro-crack if two condi- 
tions are met: (i) the local stress reaches cry and (ii) the 
SS is located at a distance from the crack tip larger than 
da 01- The nucleation of a micro-crack is assumed to 
be accompanied by an excluded volume V where stress is 
screened i.e. no SS can acivate anymore. In the follow- 
ing, <jy, e, da and V are taken as constants through- 
out the material. Three cases should be considered: 

(I) - At the onset of crack propagation, all the volume 
within Rc{v = 0) = K^jaa^ contributes to the fracture 
energy V{v = 0) = KHE. 

(II) - For < '^a, no micro-crack nucleates and Rc{v) = 
Kd{c^ v)'^ /aay < da- The dynamic SIF is then Kd{c, v) = 
k{v)K{c) 111] where k{v) 2^ (1 — v j cr) j ^/Y^^vJcb is 
universal and is the dilatational wave speed (here 
CD=2010±60 m.s-^). The volume scanned by the FPZ 
when the crack surface increases by S is 2Rc{v)S. The 
dissipated energy T{v)S is given by -\-2eRc{v)S where 
7 is the Griffith surface energy. Since T{v = 0) = K^/E^ 
one finally gets for v <Va'' 



T{v) 



a- 



E 



(1 



^2 



with a = 



2eE 



acr 



(2) 



This predicted linear dependence of F with /E for 
V < Va''^^ in agreement with measurements (Fig. [21 inset). 
A linear fit to the data (correlation coefficient i?=0.985) 
gives a=1.17±0.05 and / E={).?>±{).2 kJ/m^, where ± 
stands for 95% confidence interval. The latter value 
is compatible with the measurements of the fracture 
energy at crack initiation. By combining Eqs. ([!] ) 
and (|2j), one gets a prediction for the T{v) curve [22| 
that reproduces very well the low velo city regime in 
Fig. [21 Extrapolation of this regime [22| exhibits a 
divergence of the dissipated energy for a finite veloc- 
ity v'^={a - 1)crCd / [acD - cr) :^200 m.s"^ c^^.2?>cr, 
slightly larger than Va- In the absence of micro-cracks, 
this velocity v'^ would have therefore set the limiting 
macroscopic crack velocity. 

(Ill) - For V > Va^ Rc{v) > da i.e. micro-cracks start to 
nucleate. The surface density of micro-cracks p{v) is then 
equal to the number of activated SS beyond da per unit 
of fracture area, i.e. ps {2 [Rc{v) — Rc{va)] — pV} where 
the third term in the parenthesis stands for the excluded 
sites around micro-cracks. This yields: 

Ka{vY-K^^ . 2E 



p{v) = f3- 



E 



with /3 



aa^ 1 + psV 



(3) 
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where Ka=Kd{va). This Unear relationship is in good 
agreement with the measurements for p{K'^/E) before 
the micro-branching onset, beyond which p saturates 
(Fig. HI inset). A fit to the data (i?=0.877) between 
Ka and Ki, — Kd{vb) gives P =33±3 J~^. In the micro- 
cracking regime, the local dynamic SIF Kd is not equal 
to the macroscopic one anymore, but corresponds to that 
at the individual micro-crack tips, at which the limiting 
velocity is expected to be v^^ > Va- It is then natural to 
assume that all micro-cracks propagate at the same ve- 
locity Va, which yields Kd{v) = k{va)K {2^. The energy 
T{v)S dissipated when the crack surface increases by S 
is 7^ + e [2Rc{v)S - p{v)SV], yielding: 



r(^) 



x- 



E 



with X 



1 + 



(4) 



where = Viva)- Eq. (|4|) predicts a linear dependence 
of F with K\lE^ in agreement with the measurements for 
KllE > Kl/E (Fig. El inset). A linear fit to the data 
between Ka and Ki, = gives x = 0.67 ±0.01. The 



corresponding predicted r{v) curve [22| reproduces very 
well the intermediate velocity regime Va < v < (Fig. 
[2]) and exhibits a divergence of the dissipated energy for 
^oo = cr{1 - x^ivaf) ^ 450 m.s-i {).b2cR. This 
limiting velocity is very close to the observed maximum 
crack speed in brittle amorphous materials. 

This simple scenario allows to illustrate how material 
defects control the dynamic fracture of amorphous solids 
before the onset of micro-branching. For v < the 
mechanical energy released at the crack tip is dissipated 
into both a constant surfacic energy and a volumic en- 
ergy within the FPZ, the size of which increases with 
crack speed. With this mechanism alone, the crack speed 
would be limited to a value slightly larger than Va- But 
damage spreading through micro-cracking makes possi- 
ble to observe much larger velocities: The crack propa- 
gates through the nucleation, growth and coalescence of 
micro-cracks, with a macroscopic effective velocity that 
can be much larger than the local velocity of each micro- 



crack tip |17lJ24l|. We suggest that micro-cracks in them- 
selves do not increase dissipation, but rather decrease it 
by locally screening the stress. At velocities larger than 
^5, micro-branches contribute to the dissipated energy 
proportionally to their surface [25]. We emphasize that 
the nominally brittle to quasi-brittle transition occurring 
at Va is very likely to be generic for amorphous solids 
and should therefore be taken into account in future con- 
ceptual and mathematical descriptions of dynamic frac- 
ture. In this respect. Continuum Damage Mechanics 
(CDM) [2^ initially derived for "real" quasi-brittle ma- 
terials like ceramics or concrete may be relevant to de- 
scribe fast crack growth in nominally brittle materials. In 
particular, a better understanding of the relationship be- 
tween the dynamics of propagation of both the individual 
micro-cracks and the macroscopic crack is still needed. 
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